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o: 

<N ; ABSTRACT 

a- 

We report a result of an interferometric observation toward an SiO maser 

^ ; source, IRAS 19312+1950, in the HCO+ J = 1-0 line with the Berkeley-IUinois- 

Maryland- Association (BIMA) Millimeter Array. In the spatially integrated spec- 

>- ■ trum of HCO^, two kinematic components were seen: a strong line with a narrow 

width (~ 2 km s~^, narrow component), and a weak symmetric line with a broad 

CN ' width (~ 60 km s~^, broad component). The line profile of HCO"'" are remi- 

(^ ■ niscent of that of CO. In the integrated intensity map, we found a pronounced 

^ ■ bipolar shape consisting of lower- velocity parts of the broad component. The 

■^ ' higher- velocity part of the broad component originated in a relatively small re- 

Q^' gion (< 3"). The spatial structure of the narrow component clearly correlated 

Q ■ with near-infrared structure. The position-velocity diagrams indicated presence 

-<— > ' of a bipolar outflow with an expansion velocity of about 10 km s" . On the basis 

ci ■ of the present results, we suggest that the nature of the bipolar flow seen in IRAS 

K^ ■ 19312+1950 is explained by hydrodynamical interaction between an AGE wind 

^ ■ and ambient material with axial symmetric structure. 
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Introduction 



IRAS 19312+1950 is an SiO maser source discovered by Nakashima & Deguchi (2000) 
in their SiO maser surveys toward selected IRAS sources in the Galactic plane (Nakashima 
& Deguchi 2003a,b). Nakashima & Deguchi (2000) first suggested that IRAS 19312+1950 
is an AGB/post-AGB star evolved from a relatively massive progenitor on the basis of its 
prominent bipolar shape seen in near-infrared images and the low color temperature of the 
dust envelope. In fact, SiO maser sources are usually identified as a late-type star with active 
mass-loss. Subsequently, they have found two remarkable characteristics of this object: (1) 
two kinematic components in molecular line profiles (strong line with a narrow width of < 
2 km s^^; weak-symmetric line with a broad width of > 60 km s~^, Nakashima et al. 2004) 
; (2) a rich set of molecular species, (Nakashima et al. 2004; Deguchi & Nakashima 2003). 
In addition, recent near-infrared observations (Murakawa 2004, in private communication) 
suggested that there is a ring-like structure with a size of about 10" around the central star. 
Some of the characteristics found in IRAS 19312+1950 are often seen in dark clouds (or 
YSOs). However, secure examples of YSOs emitting SiO masers are quite limited; actually, 
only three: Ori IRc 2, W51 IRs 2 and Sgr B2 MD5 (cf., Hasegawa et al. 1986). These YSOs 
with SiO masers lie in the extreme star forming regions (in giant molecular clouds), which 
are clearly identified as largely extended nebulae in infrared images, whereas no clear star 
forming activity is seen around IRAS 19312+1950 in infrared images. Thus the evolutionary 
status of IRAS 19312+1950 is not yet definitely known. 

To investigate the nature of this highly unusual SiO maser source, we are currently con- 
ducting high-spatial-resolution interferometric observations using the BIMA array in several 
molecular rotational lines detected by our single-dish observations (Nakashima & Deguchi 
2000; Nakashima et al. 2004; Deguchi & Nakashima 2003). In this paper, we report the first 
result of the interferometeric observations in the IICO+ J = 1-0 line, and the detection of a 
bipolar fiow, which might be interpreted by interaction between an AGB wind and ambient 
material. 



2. Observation and Results 

Interferometric observations of IRAS 19312+1950 were made with the Berkeley-Illinois- 
Maryland- Association (BIMA) Millimeter Array from December 2003 to January 2004. The 
instrument is described in detail by Welch et al. (1996). We observed the HCO+ J = 1-0 
line at 89.188526 GHz with the BIMA array consisting of 10 elements in one configuration 
(B-array). The observations were interleaved every 25 minutes with the nearby point source, 
1925+211, to track the phase variations over time. The absolute fiux calibration was de- 
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termined from observations of Uranus, and is accurate to within 20%. The final map has 
an accumulated on source observing time of about 15 hrs. Typical single sideband system 
temperatures ranged from 200 K to 300 K. The velocity coverage was 380 km s~^, using three 
different correlator windows with a bandwidth of 50 MHz each. The velocity resolution was 
1.3 km s-^ The phase center of the map was 19'*33™24.4^ 19°56'54.8" (J2000) correspond- 
ing to the IRAS position of this object. Data reduction was performed with the MIRIAD 
software package (Sault et al. 1995). Standard data reduction, calibration, imaging, and 
deconvolution procedures were followed. Robust weighting of the visibility data gave 3.7" 
X 2.5" CLEAN beam with a position angle of 23.0°. The r.m.s. noise per 1.0 km s^^ is 
3.8x10^2 jy beam-^ 

Figure 1 shows the spectrum of the HCO+ J = 1-0 line. A remarkable feature seen 
in the spectrum is a strong line with a narrow width (~ 2 km s~^) peaked at Visr ~ 38 
km s^^. A weak symmetric line with a broad width (~ 60 km s^^) is also seen in the 
velocity range of Visr ~ —10 — 70 km s^^. For convenience, we call the former the "narrow 
component", and the latter the "broad component" in this paper. The narrow component is 
roughly centered on the broad component. Absorption features are seen beside the narrow 
component, suggesting that cold gas components lie in the foreground. The line profile of 
HCO"'" is reminiscent of that of the CO J = 1-0 (and J = 2-1) line (cf., Nakashima et 
al. 2004). Similar line profiles including the two kinematic components (narrow and broad 
components) have been reported in some AGB stars (ex., EP Aqr, X Her and RV Boo, 
Kahane & Jura 1996; Kerschbaum & Olofsson 1999; Knapp et al. 1998). Kahane & Jura 
(1996) reported on a bipolar flow consisting of the broad component seen in X Her, and 
Knapp et al. (1998) suggested a scenario where episodic mass-loss plays a role to explain the 
kinematic components seen in these AGB stars. Figure 2 shows velocity integrated intensity 
maps in three different velocity ranges corresponding to the narrow component (broken 
contour) and the blue- and red-shifted wings of the broad component (thin and thick solid 
contours, respectively). The map is superimposed on a combined near-infrared image (J, 
H, and /C-bands) taken by the CIAO camera on the SUBARU telescope (courtesy of Koji 
Murakawa, Motohide Tamura and the CIAO group in the National Astronomical Observatory 
in Japan). The velocity ranges used for computing integrated intensities are indicated in 
Figure 1 and the caption of Figure 2. The source is clearly resolved by the synthesized beam 
(shown at the lower-right in Figure 2). Though we observed with only one configuration, 
most part of the flux emitted from the source is thought to be detected according to the 
upper limit of the source size (~ 30") determined by our single-dish observation (Nakashima 
et al. 2004). As the most pronounced feature in the map, we can see a bipolar shape 
consisting of the thin and thick contours. The apparent axis of the bipolar shape is close to 
the NNE — SSW direction. If the bipolar shape originates in a bipolar outflow, the lower 
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side (southern component) should be closer to us, and the upper side (northern component) 
should be further from us. Interestingly, the structure of the blue-shifted wing (thin contour) 
of the broad component shows two extended tails: the strong tail to the southeast, and the 
weak tail to the east exhibiting a " lying- Y" shape. The structure of the red-shifted wing 
(thick contours) extends only to the northeast. On the contrary, the structure of the narrow 
component (broken contour) is somewhat complicated. Three intensity peaks are seen in the 
broken contour at the points of 3", 4" and 8" from the map center to the northeast, northwest 
and east, respectively. Although the positions of these peaks of the narrow component do 
not correspond to bright regions seen in the near-infrared image (background of Figure 2), 
it looks as though the global structure of the narrow component clearly correlates with the 
near-infrared structure. From the apparent size of the HCO^ structure and the assumption 
of the density, we can crudely estimate the mass of the source. If we assume that the source 
is a homogeneous sphere with a diameter of 10", the mass is estimated to be 4-31 Mq, where 
we used the distance (2.5 - 5.5 kpc) estimated by Nakashima et al. (2004) and the assumed 
density of 10^ cm~^ traced by the HCO"'" line. 

To investigate the kinematic structure, we made position-velocity (p-v) diagrams at 
various cuts. Figure 3 shows the p-v diagrams at selected cuts. In the top four diagrams 
in Figure 3, velocity channels are averaged over 2 km s~^ intervals. The cuts used for the 
diagrams are indicated in Figure 2 as the broken arrows. The directions of the arrows in 
Figure 2 represent the positive direction of the offset axes in Figure 3. The origins of the 
cuts A, B, C and D are taken at the phase center. The origin of the cut E is indicated by the 
filled dot in Figure 2. The cut "A" corresponds to the apparent axis of the bipolar shape. 
The cut "B" represents the perpendicular cut to the cut A. The cut "C" corresponds to the 
major axis of the ring-like structure seen in near- infrared images. The cut "D" represents 
the perpendicular cut to the cut C. The cut "E" is mentioned later. Although the global 
structure seen in Figure 3 is highly clumpy, the spatial size (vertical size of structure) of 
the source tends to be smaller in higher-velocity regions. Intensity peaks lie nearly on the 
origins in vertical axes in velocity ranges, Visr ^ 20 km s~^ and Visr ^ 50 km s"^. In addition, 
maximum velocity widths (represented by the maxim widths of the lowest contours in the 
horizontal coordinates) do not exhibit any noticeable variations with respect to the position 
angle. These facts mean that the higher velocity part {\Visr — Kj/s| ^ 15 km s~^; here, we 
assume that the system velocity is Vsys ~ 37 km s^^) of the broad component is not resolved 
by our synthesized beam, and should originate in a relatively small region (< 3"). On the 
other hand, variations of structure are found in lower velocity ranges {\Visr — Vsys\ ^ 15 km 
s~^) of the broad component especially in the blue-shifted side (20 km s^^ < Visr ^ 35 km 
s~^). In fact, intensity peaks at Visr ~ 25 — 27 km s^^ clearly shift to the negative direction 
in the offset axes. These peaks correspond to the southern and northern tails of the blue- 
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shifted wing (cf., thin contour in Figure 2). The 27 km s~^ peaks seen in the cuts A, C and 
D correspond to the southern tail, and the 25 km s~^ peak seen in the cut B corresponds to 
the northern tail. On the contrary, in the red-shifted side, no clear variations of structure 
are seen except for an extension to the positive direction in the offset axis in the cut A (and 
possibly in C) at Visr ~ 47 km s~^. We also made p-v diagrams without channel binding 
to check the motion of the narrow component in several cuts. Finally, we found the most 
pronounced variation of structure in the cut "E" corresponding to the tail growing to the 
southeast from the intensity peak at (DEC, RA)oflfsct =(1-1") —2.5"). The p-v diagram in 
the cut E is shown in the bottom panel in Figure 3, and there is clear variation of structure 
as a function of radial velocity. The peak has another tail extended to the southwest. The 
southwest tail also exhibits a systematic variation of structure in a p-v diagram, though 
the diagram is not shown in this paper. In any other cuts, no clear variation of velocity 
structure is found in the narrow component. We are likely to see complicated variation of 
the narrow component in the top four panels in Figure 3. These structures are explained 
well by a superposition of the two flows mentioned above. 



3. Discussion 

The bipolar flow seen in IRAS 19312+1950 consists of the lower- velocity part of the 
broad component, whereas the higher- velocity part originates in a small region with a size 
of less than about 3". Usually, the expansion velocity of bipolar flows seen in late- type 
stars (especially in post-AGB stars) tends to increase with a distance from the central star 
(frequently called "Hubble type flow" in this area, cf., Balick & Frank 2002), the nature 
of the bipolar flow seen in IRAS 19312+1950 is in disagreement with this point. In our 
opinion, the characteristics of the bipolar flow are explained by hydrodynamical interaction 
between a spherical outflow lying at the center and ambient material with axial symmetric 
structure if the spherical outflow expelled from the central star is distorted to a bipolar 
shape by the interaction. In such a case, the slow expanding velocity of the bipolar flow 
seems to be somewhat strange, but the slow velocity would be explained by a projection 
effect if the flow has an inclination angle of 20° (where we assume the velocity of the bipolar 
flow equals to a half of the maximum width of the broad component, i.e. ~ 30 km s^^). 
Interestingly, in Figure 2, the bipolar shape precisely escapes from dense regions of the narrow 
component, and the red-shifted component of the bipolar flow lies in the interspace between 
two remarkable intensity peaks of the narrow component. In addition, the absorption feature 
beside the narrow component seen in Figure 1 might be explained as the near side of the 
ambient material (with axial symmetric structure). 
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If the higher- velocity part of the broad component is a spherical outflow from the central 
star, the flow should be identified as an AGB outflow, because the spherical outflow is the 
typical nature of AGB envelopes, and also because the expanding velocity of Vexp ~ 30 km 
s~^ is very reasonable as that seen in the super-wind phase at the late-AGB stage. Our recent 
CO data taken by the BIMA array also support the spherical property of the inner part of 
the nebulosity. However, a problem is the somewhat strong HCO"'" intensity (especially in 
the outermost parts of the envelope) if the central star is an 0-rich AGB star, because the 
HCO"*" line is usually weak or null in AGB envelopes except a few cases (Nguyen-Q-Rieu et 
al. 1988; Deguchi et al. 1990; Cox et al. 1992). From the viewpoints of theory, the HCO+ 
molecule can be produced by photo-chemical process in an 0-rich AGB envelope up to the 
radius ~ 2 x 10^^ cm (e.g., Mamon, et al. 1987, ; Mamon et al. took account of ionization by 
galactic and stellar UV photons, and high-energy particles in their calculation.). However, 
we still need additional chemical process to explain the effective formation of HCO"*" in the 
outermost parts of the envelope, because it is most likely that the outermost part of HCO"*" 
emission is further than 10^^ cm from the central star if the central star is an AGB star 
(estimated distance: 2.5-5.1 kpc, see Nakashima et al. 2004). An immediate possibility to 
explain the highly extended HCO"*" structure would be "dissociative shock". Actually, most 
AGB envelopes exhibiting the HCO"*" emission show signs of shocks (cf.. Cox et al. 1992). 
The HCO"'" emission in OH 231.8-1-4.2, which is a well-known example of an SiO maser 
source with a bipolar shape, is also explained by shocks (Sanchez Contreras et al. 2000). 
The shock ionization also supports our idea to explain the kinematic structure of the bipolar 
flow, because shock fronts are likely to be produced between an AGB wind and ambient 
material. 

Remaining problems in the interpretation as an AGB star are excess of mass and mass- 
loss rate. Although we estimated mass to be 4-31 Mq in a previous section, the mass will be 
over 8 Mq if the distance exceeds 3 kpc. In such a case, the nebulosity cannot be explained as 
material expelled from a central AGB star. If we assume a spherical flow and density of 10^ 
cc^^ at 10^'^ cm from the central star, the mass-loss rate is estimated to be about 10"'^ M© 
yr^^ in terms of the expanding velocity of the broad component (~ 30 km s^^). This mass-loss 
rate is slightly larger than a typical value of AGB stars (~ 10""^ Mq yr^^). The keys to check 
the AGB possibility would be near-star kinematics and isotope ratio (ex., ^^C/^^C), because 
AGB stars and YSOs are known to be have significantly different kinematics especially in 
a near-star region, and because ^^C should be enhanced in a circumstellar envelope if the 
central star lies at the late-AGB stage. 
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4. Summary 

In this paper, we reported a result of an interferometeric observation of an SiO maser 
source, IRAS 19312+1950, in the HCO^ J = 1-0 line. We found two kinematical components 
in the spectrum including a strong-narrow line (narrow component) and a weak-broad line 
(broad component). Velocity integrated intensity maps showed a pronounced bipolar shape 
consisting of lower-velocity parts of the broad component. The higher- velocity parts of the 
broad component were not resolved by our synthesized beam, and should originate in a 
relative small region with the size of less than ~ 3". The spatial structure of the narrow 
component clearly correlated with the near-infrared structure. In p-v diagrams, a systematic 
variation of structure was found in the lower- velocity part of the broad component, indicating 
a presence of a bipolar outflow. We suggest that the characteristics of the bipolar flow seen in 
IRAS 19312-1-1950 are explained by interaction between an AGB wind and ambient material 
with axial symmetric structure. 

The authors are grateful to Kee-Tae Kim, James R. Forster, Hiroshi Imai, and Hideyuki 
Izumiura for stimulating discussion, and to Anthony Remijan for help with the BIMA ob- 
servations and data reduction. We would like to also thank Douglas Friedel for help in 
correction of the English, and the anonymous referee for useful comments. This research 
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Fig. 1.— Spatially integrated spectrum of IRAS 19312+1950 in the HCO+ J = 1-0 line. 
The integrated area is a circle with a diameter of 15". The vertical solid lines represent the 
velocity ranges of the blue- and red-shifted wings of the broad component used for Figure 2. 
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Fig. 2. — HCO^ integrated intensity map (contours) superimposed on a near-infrared com- 
posite (J, H and K bands) image (grayscale) taken by the SUBARU telescope. The ranges 
of velocity integration for the thick, thin and broken contour maps are 0-26 km s~^, 47 
- 60 km s~^ and 37 - 38 km s~^, respectively. The contour levels start from a 5 a level 
, and increase by every 0.5 a for the broad component and by every 1 a for the narrow 
component. The 1 a levels of the thin, thick and broken contour maps are 7.3 x 10~^ Jy 
beam~^, 1.0 x 10~^ Jy beam~^ and 2.7 x 10~^ Jy beam~^, respectively. The synthesized 
beam is indicated in the lower right corner. The broken arrows (A, B, C, D and E) denote 
"cuts" used for p-v diagrams in Figure 3. The filled dot represents the origin of the cut E. 
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Fig. 3. — Position-velocity diagrams along the cuts indicated in Figure 2. The velocity 
channels are averaged over 2 km s~^ intervals except for the bottom panel. In the top four 
diagrams, the contours are drawn from 40 mJy beam~^ with increments of 10 mJy beam~^ 
between 40 and 80 mJy beani"^, and 20 mJy beam~^ between 80 and 160 mJy beani"^. The 
bottom panel is the p-v diagram along the cut "E" without velocity channel binding. In the 
bottom panel, the contours start from 60 mJy beam~^ with increments of 20 mJy beam~^. 



